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Lattice QED and Universality of the Axial Anomaly 
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We give a perturbative proof that U(l) lattice gauge theories generate the axial anomaly in the continuum limit 
under very general conditions on the lattice Dirac operator. These conditions are locality, gauge covariance and 
the absense of species doubling. They hold for Wilson fermions as well as for realizations of the Dirac operator 
that satisfy the Ginsparg- Wilson relation. The proof is based on the lattice power counting theorem. The results 
generalize to non-abelian gauge theories. 



1. Introduction 

In the continuum formulation of QED (or 
QCD) it is well known that all known gauge in- 
variant regularization schemes yield the ABJ ax- 
ial anomaly. A non-perturbative regularization 
is provided by the lattice. Any candidate for a 
lattice discretization of the QED action should 
reproduce the axial anomaly in the continuum 
limit. A naive symmetric discretization of the 
continuum Dirac operator fails to generate the 
anomaly, but is also plagued with the fermion 
doubling problem. The doubling problem can be 
eliminated by introducing an explicit, naively ir- 
relevant, chiral symmetry breaking Wilson term 
in the action. This circumvents the Nielsen Nino- 
myia theorem [[lj and leads to the correct anomaly 
in the continuum limit 0,0]. 

The question arises whether this is a particu- 
lar property of the action chosen, or whether the 
anomaly will be correctly reproduced by any lat- 
tice regularization of the action satisfying some 
very general conditions. In fact, as we will show, 
any lattice action which is i) gauge invariant, 
ii) possesses the correct continuum limit, iii) is 
local in some generalized sense, and iv) is free 
of doublers, will necessarily reproduce the axial 
anomaly in the continuum limit. An important 
ingredient which will be needed for the proof is 
the general form of the lattice axial vector Ward 
identity. Its precise structure, which will depend 
on the particular discretization of the action cho- 



sen, need not be known, but only very general 
properties thereof. Our general proof will include 
the case of Ginsparg- Wilson fermions. 

2. Ward Identity 

Consider the following general form for the 
fermionic contribution to the lattice action 



ip(x)(Du(x,y) +m)ip(y) . (1) 
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Here U stands for the collection of link variables, 
Du is the Dirac operator, and J2 X = S n fl4 ' 
where n label the lattice sites. Du can be decom- 
posed into a part anticommuting with 75, which 
in the continuum limit takes the form 7 • d, and 
a term commuting with 75, which vanishes in the 
continuum limit. Hence the general structure of 
the lattice axial vector Ward identity in an exter- 
nal gauge field will be of the form 

(d;jf t (x)) u = 2m{j 5 (x)) u + (A(x)) u (2) 

where <9*(<9 M ) stands for the dimensioned left 
(right) lattice derivative. j 5 (x) — ip(x)^i^{x), 
and jn( x ) nas the correct continuum limit. A (a;) 
is an "irrelevant" operator which vanishes in the 
naive continuum limit. For Wilson fermions jfj 
and A have been given in [pj. Furthermore, it was 
shown in that the anomaly is generated in 

the continuum limit by the naively "irrelevant" 
contribution < A(x) > . For Ginsparg- Wilson 
fermions the Dirac operator in (Q) satisfies the 
Ginsparg Wilson relation PJ 

75-D + D75 = aD^ 5 D , 
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where a is the lattice spacing. The action pos- 
sesses an exact global axial symmetry jj| for 
m = under the transformations 

*#c)=iw7ji[(l-|l>M(aO (3) 

<^(x)=iu#(l-!l))](z)7 5 • (4) 

According to the lattice Poincare Lemma Q there 
exists an axial vector current associated with this 
symmetry. Under a local transformation with 
lj replaced by u){x), the variation of the con- 
tribution to the fcrmionic action involving the 
Dirac operator in ([!]) can be written in the form 
5Sd = —'iJ2x UJ ( x )^^fi( x )- Taking into account 
the variation of the fermionic measure one then 
finds that the axial Ward identity is of the form 
, where the irrelevant A operator is given by 

A GW (x) = a$D)(x)K(Dip)(x) . 

Consider now the Ward identity (^) in momentum 
space. Define < j^{q) >u by 

< 3%{x) >u= f\ ^jL < #( g ) >v e i( 9 -*+-A/2) 

Expanding U^x) — exp[igaA^{x)] in the poten- 
tials, < i^(q) >u will have the form 

x^i/i—un (fcii ■ ■ ■) k n )A ! ; 1 (fei) • • ■ A Vn (k n ) , 

with analogous expressions for < j 5 (q) >u and 

< A(g) >u- The Ward identity in momentum 
space then reads 

^r^ } ({fc l }) = 2mr^ s} ({fc l }) + r^ ) } ({fcj) ,(5) 

where q^ = - sin ^ . 

3. Theorem 

We now state the central theorem of this paper 
and then give the proof: 

Any lattice discretization of the QED action 
with the following properties: a) S has the correct 
continuum limit; b) S is gauge invariant; c) The 
Dirac Operator is local and d) The free Dirac op- 
erator ZH°)(p) is free of doublers, reproduces the 
axial anomaly in the continuum limit. 



As we shall see below, the explicit form of the 
axial vector current and of the irrelevant operator 
A(x) is not required. Let us consider a) to c) in 

turn: 

1) S has the correct continuum limit: This en- 
sures that j~ posesses the correct naive contin- 
uum limit. 

2) Gauge invariance: Gauge invariance tells us 
that if 0{ip, ip,A) is a gauge invariant operator, 
then its external field expectation value satisfies 

< OU>,$,A u ) > A ^< 0(ij,^,A) > A , 

where A"(x) — A a (x) + d^ui(x). It follows that 

or in momentum space 

Y / CkU^ ) -,Aku---,k n ) = . (6) 

3) Locality of the Dirac Operator. The Dirac op- 
erator Du(x,y) can be formally expanded in a 
power series in the gauge potentials: 

Du(x,y) = En, W)Zj ^. D )^... tin {x,y\z 1 ■ ■ ■ z n ) 
xA w (2i) ■ • • A^ n (z n ) . 

The coefficient functions D^".. tJ[/n are assumed to 
vanish exponentially fast with the separation be- 
tween any pair of lattice sites with a decay con- 
stant of the order of the lattice spacing (this 
also holds for Ginsparg Wilson fermions for suf- 
ficiently smooth gauge field configurations 0. 
As a consequence T l %... Vn ({k i }) , ^...^({fej}), 

and T^).. Vn ({ki}) can be Taylor expanded around 
zero momenta. It then follows from (|^) that 

%Tvi,-,v n {k\r ■ "jfcn) = , (7) 

where T n denote the Taylor expansion up to order 
n in the momenta 

4) The free propagator (p) is free of doublers. 
As a consequence the lattice power counting rules 
H and the Reisz theorem || applies. This theo- 
rem allows one to take the naive continuum limit 
of a lattice integral, if all the lattice degree of di- 
vergences of the integrand, including the measure, 
are negative. 
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Proof of the Theorem 

In the following we only consider the axial vec- 
tor Ward identity for two external photons, since 
it is the only one which is anomalous. It has the 
form 

q^M*}) = 2mr^^({Ai})+rW ,(8) 

where i = 1,2. By power counting degrT 5fJ " = 
degrT 5 — 1, while the irrelevant contribution has 
degr r( A ) < 2. Here degrF denotes the lattice 
divergence degree of F ||. 

We first show that all the vertex functions pos- 
sess a continuum limit. Consider e.g. T 5fi , and 
write it in the form 



rfA = (1 - TOrf/ 1 ,, + TxTl» 



(9) 



The first term on the RHS has an degr < 0. 
Hence making use of the Reisz theorem, its con- 
tinuum limit is given by the Taylor subracted con- 
tinuum one loop integral with a 7^75 insertion. 
This is just the usual continuum triangle graph, 
Taylor subtracted to first order in the momenta. 
The second term on the RHS of (^j) vanishes be- 
cause of gauge invariance (^) . The a — > limit of 
„ 2 can be calculated in a similar way. Hence 
we immediately arrive at the result that the con- 
tinuum limit of Tif} 2 in (||) is universal, and given 
in terms of convergent continuum Feynman inte- 
grals. This limit can be easily computed. To this 
effect decompose P^f^ as follows: 

vi% = (1 - r 2 )rifj 2 + 7ir(f) 2 + (r 2 - r^rw 

where T 2 denotes Taylor subtraction to second 
order in the momenta. Because of gauge invari- 
ance (fjj) the second term on the RHS vanishes. 
According to the Reisz theorem M the first term 
vanishes for a — > since degr[r( A )] < 2, and 
vanishes in the naive continuum limit. Hence we 
conclude that 

lim Y^l (h ,k 2 ) = km (T 2 - 7i)4 A » (h , k 2 ) 

We now make use of the Ward identity (||) to 
obtain a simple expression for the RHS. Thus ap- 
plying the T 2 — Ti operation to both sides of (||) , 
making use of gauge invariance ((7J) and of the 
Reisz theorem, one readily finds that 

dH 



where I^ lVa (i, ki,k 2 ) is the integrand of the con- 
tinuum Feynman integral for the triangle graph 
with a 75 insertion. Hence the anomalous contri- 
bution reads: 



-16e e^ u ^pk\ a k 

1 

2^2 



d A l 



l ' ,n,2p I (2tt) 4 (^ 2 +m 2 )3 
~z j ^ , ^ fiu pak\pk 2a 



p,a 

We have thus shown that, irrespective of the de- 
tailed structure of the lattice current and of the 
irrelevant operator A in the Ward identity (||), 
the correct anomaly is generated by the irrelevant 
A-contribution in the continuum limit. This limit 
is universal and given by the second order Taylor 
term in the expansion of the continuum triangle 
graph with a 75 vertex insertion. This expression 
has precisely the same form as that obtained in 
the continuum formulation using the Pauli-Villars 
regularization scheme, except that here the Pauli 
Villars mass is replaced by the fermion mass. The 
integral is however independent of the mass. Fur- 
thermore we have seen that all terms in the Ward 
identity (||) posssess a continuum limit. The limit 
of the first two terms is however not given by the 
naive limit of the triangle graphs, but actually by 
their BPHZ subtracted form. 
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